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Q\ • A young open cluster is a 2-phase system. 
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We present an approach which consists in a (simplified) model where stars (N-bodies) 
move within a gaseous spherical molecular cloud. The two components influence each 



• an ensemble of stars move in a gaseous medium (the mother molecular cloud). 

The dynamics and thermodynamics of the system, and so its 
-evolution and final fate (is it stable or unstable?) 
strongly depends on the mutual feedback between gas and stars. 



other through 

- gravity and mass loss. 

Among other results (role of IMF, SFE, stellar background, etc., see Conclusions), 
we find that a significant fraction of small clusters can be destroyed even 

- before SN explosion, 
when a significant amount of massive stars are present. 

THE MODEL 

After the Lada, Margulis and Dearborn (1984, ApJ 285, 141 LMD) work not 
much has been done to study quantitatively the early evolution and fate of stellar 
clusters embedded in their mother cloud, following numerically the N-body dynamics 
of stars moving in a (dispersing) gaseous cloud. The LMD model was not fully self- 
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consistent, for the gas was assumed to expand with an assumed time law; by the way, 
this work gave relevant information on the capability of a stellar system to remain 
bound after gas removal in dependence on the star formation efficiency (SFE). 

- An answer to the crucial question: 

• what conditions on IMF and on SFE allow a small cluster, emerging from a molec- 
ular cloud, to remain bound? 

- necessarily implies that the mutual feedback between gas and stars is taken 
into account. 

To get really reliable results one should couple an N-body code to a fully hydro- 
code to model the radiative and mechanical interaction between the stellar and 
gaseous phases. This has been partially done (see Capuzzo-Dolcetta and Di Lisio, 
SPH in Astrophysics, 1994, Mem.S.A.It., 65, 1107), and is the target of future work 
(Capuzzo-Dolcetta, Di Lisio, Navarrini, Palla, in prepaparation) . 

Results good to order of magnitude, can however be obtained with the present 
model, which treats the coupled dynamics and thermodynamics of stars and gas in 
a cluster with the following: 

approximations 

• the gas cloud evolves in time keeping spherical shape and a spatially uniform density. 

• the gravitational force exerted by stars on the cloud is approximated. 



The Equations 



The relevant equations are: 



a,- = 



drrij Vi 



i = l,...,N 



mi dt m 

R = Rg + R p + -R* + R m l + Rvr 
KU = U p + U vr + Umi + U SN 

which is a 6N + 3 order system submitted to the appropriate initial conditions. 

- di is the i-th star's acceleration, 

- R is the gas-sphere radius, 

- U is the gas internal energy. 
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In the gas motion equation: 

« GM , ,„ 
9 = — W ^ ~~ gravity) 

= ^ ^Zr (pressure field) 

MgUg 

1 * 

= ^ fi (stellar - gravity) 
i=i 

R m i = 2_. ( 1 ~~ 2^ J -Rmii (stellar mass — loss) 

i= i V / 



M + M: 



{ ■ N \ 

- I Mi? + ^^miri J (vzoZ. relax.) 

* V *=i / 



• /j is an approximation of the force exerted by the i-th star on the gas cloud. 

• v m i i is the z-th star wind speed (taken from the literature). 

In the gas energy equation: 

rr 9 / r 

Up =-7(7 — l) - ^ - (pressure heating) 

3 M ■ ( ■ N \ 
U vr =-k vr — 7~rR I MR + ^^mifi J (viol, relax.) 



5 U 'M + M* 

N 



i=l 



Usn =S(t — tsN)esN (SN contribution) 



(star mass — loss) 
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Parameters of the models 

Stars are initially uniformly distributed in space and velocity in a sphere of radius 
R*o with velocities to satisfy the given virial ratio (here assumed =1). 

The relevant initial parameters are: 

Stars: 

N = number of stars 
R*o = initial cluster radius = 1 pc 
IMF oc m" Q , 0.2 < m/M© < 20 
local SFE = e 

chemical composition = (X, Y, Z) — (0.7, 0.27, 0.025) 

2 x kinetic energies 

virial ratio = uq = ; = 1 

potential energy 

Gas Clump: 

Ro = initial radius = R*o 
po = initial gas density = 500 M & /pc 3 
Ro = initial collapse velocity = 
Uq = initial internal energy 
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CONCLUSIONS 

• A small cluster (N < 100, 0.2 < m/M & < 20) embedded in a gas clump of typical 
density p ~ 500M o /pc 3 

is lost to the background 
at the time t ~ lOMyr when first SN explode (see Fig. 1), if SFE<0.4 . 

• When SFE> 0.4, the cluster resists to the explosion whenever the IMF is not 
biased towards large masses. 

• When the exponent of the IMFoc m~ a is sufficiently negative (a < —2) and 
SFE>0.35, the gas cloud and its embedded cluster are disrupted by powerful 
stellar winds in a shorter time (few Myrs). In this case, 

just very high SFE allow the cluster to survive (Fig. 2, Fig. 3 and Fig. 4)- 

• The capability to distinguish a small cluster over a background strongly depends 

on the cut-off' magnitude (the density contrast falls of a factor 50 when V cu t is changed 
from 14 to 18 !). This means that intrinsecally bound cluster can be misconsidered 
as unbound just because they are observed over a too crowdy background: 

so any practical definition of cluster lifetime must take into account the bachground 
over which the cluster is projected (Fig. 5 and Fig. 6). 
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Figure Captions 



Fig. 1 : Super-Nova explosion time vs star mass (chemical composition X = 
0.7, Y = 0.27, Z = 0.025). 

Fig. 2 : Gas cloud dissolution time vs SFE for the IMF exponent a = —2. 

Fig. 3 : Rough delimitation of regions of bound and unbound clusters in the (a,e) 
plane, where a = IMF exponent, e =SFE. 

Values of e above the upper horizontal line have not yet been investigated. 

Fig. 4 : For the models whose a, e, N are labelled on the top : 

solid lines — cluster Lagrangian radii of 25%, 50%, 75% and 100% of the mass 

dashed lines - - gas cloud radii. 

Bottom panels are the enlargements of the upper ones. 

Fig. 5 a,b : Time evolution of the cluster density contrast Ap/p^g = (p — Pb g )/pb g - 
V cu t is the lower luminosity cut-off of the background, which is taken at latitudes 
6 = 0° and b = 45° (upper and lower curves, respectively, in each panel). Horizontal 
lines correspond to Ap/p^g = 1, below which the cluster is undistinguishable over 
the estimated background. The various cases studied label the panels. 

Fig. 6 : Enlarged view of part of Fig. 5b, to show clearly the transition from 
"visible" to "unvisible" cluster. 
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